Periodic control laws for bilinear quantum systems with discrete spectrum 
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Abstract — We provide bounds on the error between dynamics 
of an infinite dimensional bilinear Schrddinger equation and 
of its finite dimensional Galerkin approximations. Standard 
averaging methods are used on the finite dimensional approx- 
imations to obtain constructive controllability results. As an 
illustration, the methods are applied on a model of a 2D rotating 
molecule. 

I. INTRODUCTION 

A. Physical context 

The state of a quantum system evolving on a finite dimen- 
sional Riemannian manifold f2, with associated measure p,, 
is described by its wave function, that is, a point in the unit 
sphere of L 2 (il, C). A system with wave function ip is in a 

subset uj of n with the probability /|?/>| 2 d/i. 

J uj 

When submitted to an excitation by an external field (e.g. 
a laser) the time evolution of the wave function is governed 
by the bilinear Schrodinger equation 

= -\^ + V{x)^{x,t)+u{t)W{x)^{x,t), (1) 

where V, W : Cl — > R are real functions describing respec- 
tively the physical properties of the uncontrolled system and 
the external field, and u : R — > R is a real function of the 
time representing the intensity of the latter. 

B. Quantum control 

A natural question, with many practical implications, is 
whether there exists a control u that steers the quantum 
system from a given initial position to a given target. 

Considerable efforts have been made to study the control- 
lability of gj. We refer to [14], [10], [2], [9], [1], [3] and 
references therein for a description of the known theoretical 
results concerning the existence of controls steering a given 
source to a given target. As proved in [10], [8], [11], 
approximate controllability is a generic property for systems 
of the type ([TJ. 

The main difficulty in the study of ([TJ is the fact that the 
natural state space, namely L 2 (£l, C), has infinite dimension. 
To avoid difficulties when dealing with infinite dimensional 
systems, for example when studying practical computations 
or simulations, one can project system ([TJ on finite dimen- 
sional subspaces of L 2 (il, C). Obviously, a crucial issue is 
to guarantee that the finite dimensional approximations have 



dynamics close to the one of the original infinite dimensional 
system. 

C. Aim and content of the paper 

The contribution of this paper is twofold. First, in Section 
p] we provide an introduction to the class of weakly-coupled 
bilinear systems (see Definition [TJ. A feature of these sys- 
tems is that their dynamics is precisely approached by the 
dynamics of their Galerkin approximations (Proposition RJ. 
In a second part, we apply general averaging theory for the 
approximate control of finite dimensional bilinear conserva- 
tive systems using small amplitude periodic control laws. The 
method is both very selective with respect to the frequency 
(which is a good point for quantum control) and extremely 



robust with respect to the shape of the control (Section III I. 
Moreover, it provides easy and explicit estimates for the 
controllability time, the L 1 norm of the control and the 
error. Together with the results of Section |TTJ this method 
provides a complete solution for the approximate control of 
infinite dimensional bilinear quantum systems with discrete 
spectrum and time estimates. As an illustration, we consider 



the rotation of a planar dipolar molecule in Section IV 



II. WEAKLY-COUPLED BILINEAR SYSTEMS 

A. Abstract framework 

We reformulate the problem ([TJ in a more abstract frame- 
work. This will allow us to treat examples slightly more 
general than ([TJ, for instance, the example in [5, Section 
III. A]. In a separable Hilbert space H endowed with norm 
|| • || and Hilbert product (•, ), we consider the evolution 
problem 

= (A + u(t)B)1> 



d.t 



(2) 



where (A, B) satisfies the following assumption. 

Assumption 1: (A, B) is a pair of linear operators such 
that 

1) A is skew-adjoint and has purely discrete spectrum 
(— iAfe)fc e N> the sequence (\k)keN is positive non- 
decreasing and accumulates at +oo; 

2) B : H — > H is skew-adjoint and bounded. 

In the rest of our study, we denote by (4>k)kG'N an Hilbert 
basis of B such that A<f>k = —\\k<t>k for every k in N. 
We denote by D(A + uB) the domain where A + uB is 
skew -adjoint. 



Together with Kato-Rellich Theorem, the Assumption T]2 
ensures that, for every constants u in R, A+uB is essentially 
skew-adjoint on D(A) and i(A + uB) is bounded from 
below. Hence, for every initial condition ipo in H, for every 
u piecewise constant, u : t M> J2j u jX(t j ,t j+1 )(t), with 
= < < ^i < ••• < ijv+i and u ,...,um in R, one 
can define the solution 1 1-> of |2]) by 

TJVo = e (*-*i-0(*+«i-i*>o 

o e (t J -i-* 3 -2)(A+n 3 _ 2 S) Q . . . Q gtofA+noB)^^ 

for i G For a control w, in i 1 (R) we define the 

solution using the following classical continuity result. 

Proposition 1: Let u and (it n ) n gN be in L 1 (R). If for 
every t in R J u„(r)dr converges to J Q w(r)dr as n tends 
to infinity, then, for every t in R and every ip in H, 
(T""^o) n gN converges to TJ'V'o as n tends to infinity. 

B. Energy growth 

From Assumption T[T the operator iA is self-adjoint 
with positive eigenvalues. For every ip in D(A), \Aip = 
SjeN (<Aj' V')0j- F° r every s > 0, we define the linear 
operator \A\ S := (L4) s by |A|> = £ N Aj(^-, V)^, for 



E, eN Af|(^»| 2 



< 



every ip in L>(|A| S ) = {ip *E H 

+oo}. We define the s-normby \\ip\\ s = \\ \A\ s ip\\ for every?/; 
in D(|A| S ). The 1/2-norm plays an important role in physics; 
for every ip in £>(|A| 1/2 ), the quantity \(Aip,ip)\ = \\^\\y 2 
is the expected value of the energy. 

Definition 1: Let (A, B) satisfy Assumption [I] Then 
(A, B) is weakly-coupled if there exists a constant C such 
that, for every ip in D{\A\), \$t{\A\i>, Bip)\ < C\(\A\ip,ip)\. 
The coupling constant c(A, B) of system (A, B) is the 

quantity sup^^i^^o} 

The notion of weakly-coupled systems is closely related 
to the growth of the expected value of the energy. 

Proposition 2: Let [A, B) satisfy Assumption [T] and be 
weakly-coupled. Then, for every ip 6 Z>(|^4| 1//2 ), K > 0, 
T > 0, and m in i 1 ([0,oo)) for which < K, one has 

||T-(^)|| 1/2 <e c (^||Vo||i/ 2 . 

Proof: We present here a simple proof in the special 
case where ip Q belongs to D(A 2 ), u is piecewise constant, 
and D{\A\ 2 ) = D(\A + uB\ 2 ) for every u in R. This last 
equality holds for the most common physical examples. A 
general proof of Proposition [2j involving rather technical 
regularization techniques to relax this extra assumption is 
presented in [5, Appendix]. 

First note that, for every t > 0, the set D(\ A\ 2 ) = D(\ A + 
uB\ 2 ) is invariant for the unitary map ip i-> e lli - A+wBS, ip. 
Moreover, for every ip in D(\A + uB\ 2 ), the mapping t >-> 
\A + Y J uB\ 2 e t{A+uB ^ = e t( * A+uB ^\A + uB\iP is C 1 from 
[0, +00) to H, with derivative t ^ (A + uB)e t(A+uB ^\A + 
uB\ip = \A + uB\(A + uB)e^ A+uB) ip. In other words, the 
mapping t i-> e ^ A+uB ^ip is q\ from [ 0j +00 ) to D(\A + 
uB\)=D(\A\). 

Fix u : [0, +00) — > R piecewise constant, ipo in 
_D(|yl| fc+1 ) and consider the real mapping / : t M> 
(\A\ k Tfip , TJ'Vo). Since ip belongs to D(\A+u{t)B\ k+1 ), 



then / is absolutely continuous and for the argument above 
is piecewise C 1 . For almost every t, 



|(|A| fc T^ ,T^o) 



= 2U{\A\rW (h (A + u(t)B)rW ) 

= 2u(t)$i(\A\rw ,Brw ). 

Since (A, B) is weakly-coupled, one has 

\f(t)\ < 2\u(t)\\(\A\rW ,BrWo)\ 
< 2c(A,B)\u(t)\f(t). 

From Gronwall's lemma, we get 

(\A\m,m)<e 2c(A ' B ^> l{T)dT \\Mi/2 

for every ip in D(|A| 2 ). I 

C. Good Galerkin approximation 

For every N in N, we define the orthogonal projection 

7t N : ip G H i->- y ] (<j>j,i/>}<t>j G H. 



j<N 



Proposition 3: Let (A, B) satisfy Assumption [T] and be 
weakly-coupled. Then, for every n G N, N G N, {ipj)i<j< n 
in Dd^l 1 / 2 )™, and for every L 1 function u, 

c(A,B)\\u\\ L i II-/, II 

11(1(1-^)^(^)11 < 7=^' O) 

V A iv+i 

for every t > and j = 1, . . . , n. 

Proof: Fix j G {1, . . . , n}. For every N > 1, one has 



|(Id-^)T«(^)ll 



OO 

E 

rl= Ar+l 



K^n,T?(Vi))| S 



< 



A^ +1 ||TWj)llx /2 - 



< 



for every f > 0, 11X^(^)11 1/2 
2 , 2 . The conclusion then follows by 



By Proposition 

e 2c k (A,B)\\u\\ Ll 

Proposition [T] 

Remark 1: Since B is bounded, then 
\\B(Id — Wff)T^(ipj)\\ tends to as N tends to infinity 
uniformly with respect to u of L 1 -norm smaller than a 
given constant. 

Definition 2: Let N G N. The Galerkin approximation of 
Q of order is the system in H 

±=(A {N) +u{t)B {N) )x (Ejv) 

where = ttatAttat and = itnBiitn are the 

compressions of A and £? (respectively). 

We denote by X^ N ^ (t, s) the propagator of ( |S jy| > for a L 1 
function u. 

Remark 2: The operators and B^ are defined on 

the infinite dimensional space H. However, they have finite 
rank and the dynamics of (Eat) leaves invariant the A^- 
dimensional space Cm = span 1<;<Ar {0/}. Thus, (Sat) can 
be seen as a finite dimensional bilinear system in Cm- 

Proposition 4 (Good Galerkin Approximation): Let 
(A, B) satisfy Assumption [T] and be weakly-coupled. Then 



for every e > 0, K > 0, n £ N, and (ipj)i<j< n in 
D(\A\ 1/2 ) n there exists N e N such that for every L 1 
function it 

|M| L1 <X => ||T?(^-) -X^ (4,0)^11 <e, 

for every t > and j = 1, . . . , n. 

Proof: Fix j in {1, . . . , n} and consider the map t i-> 
ttjvT" (V'j) that is absolutely continuous and satisfies, for 
almost every t > 0, 



and 



(4) 



(5) 



+ U (i)^v.B(Id- 7 r w )T; i W> J ) 
Hence, by variation of constants, for every t > 0, 

7r JV T«(^-)=^(iv)(*>0)7riv^ 

+ /x ( u 7V) (t, S )7r A r J B(Id-^)T^(^)w(r)dr 
Jo 

By Proposition [5] the norm of t M> _B(Id — 
ttaOT^Oj) is less than \\B\\ X^+i^^Wih/z- Since 
is unitary, \\n N T^(^) - ^(t, 0)7^11 < 
hll^llSllA^V^^HV^II^.Then 

||Tn^-)-^)(i,0)7r^ll 
< ||(Id-^)T«(^)|| 

+\\n N rn^)~X(N ) (t,0M J \\ 



< X-^(l + K\\B\\). 



,c(A,B)K I 



This completes the proof since Xn tends to infinity as N 
goes to infinity. ■ 

III. Control of finite dimensional conservative 

BILINEAR SYSTEMS 

A. Averaging results 

In this Section, we focus on the control system (Ejv) of 
Definition [5] The matrix AW is diagonal with eigenvalues 
(— iAj)i<j<jv with Xi < A 2 < . . . < Ajy- We denote with 
(bjk)i<j,k<N the entries of B^ N \ 

Theorem 5: Let j, k be two integers such that 1 < j < 
k < N. Assume that bjk ^ and that for every I, m < N, 
\Xi - X m \ € N|Aj - Afe| implies {j, k} = {I, m} or b lm = 0. 

Then, for every u* in i 1 (R) satisfying 



Define T 



2tt 
At, — A, 



f M*(T)e i(Afc -^ )r dr ^0 and / w* (r)e i(A; - A ™ )T dT = 0, 
Jo Jo 

for every / , m such that | A j — A& | G | A/ — A m | N and bi m =^ 0, 
one has 



with 



< 



C 



|(^,X-(nT*, 0)^)| 

C 7T (joVWI 



dr 



n 2|6,-J 



Jo M*(r)e i ( A J- Afc ) T d7 

awii + iibwii 



inf a ; - 



sin ( 27r 



.a,-a„ 



ttT 



2 \bjk\ /„«*(r)e 1 (^-A*)rd-: 

Theorem [5] above states that |(</>fc, -X" 3 ^ (nT* , 0)<fij) | tends 
to one as n tends to infinity. This convergence result may be 
obtained using classical averaging theory (see for instance 
[13]). Another proof based on the particular algebraic struc- 
ture of the system is given in [6]. The estimates given here 
are a consequence of [6, Lemma 8]. 

B. Efficiency of the transfer 



Using the notations of the last paragraph, for every pe- 
riodic function u* with period T = > x we define the 
efficiency of u* with respect to the transition (j, k) as the 
real quantity: 



r o T u*(r)e i (^- A *) T dr 



/ V(r)|d7 



For every u, < E^' k \u) < 1. For every {j,k}, 
sup u E^' k ^ (u) = 1 (consider a sequence of controls that 
tends to a periodic sum of Dirac functions). An example of 
u* with zero efficiency is presented in Section IV 



An intuitive explanation of the efficiency could be the 
following: asymptotically, the L 1 -norm of the control needed 
to induce the transition between levels j and k using periodic 
controls of the form u„ is equal to -k / (2\bjk\E^' k ^ (u*)). 

C. Design of control laws 

The system |2]i being given, the design of an effective 
control law fulfilling the hypotheses of Theorem [5] is an 
important practical issue. To generate a transfer from level j 
to level k, one should chose a control u such that E^' k ^ [u) 
be as large as possible and E^ l,m '(u) be zero (or arbitrarily 
close to zero) for every I, m such that A; — X m E (Aj — Afe)Z. 

A natural choice for the control law can be t H> 
cos (\Xj — Xk\t) . The efficiency with respect to transition 
(j, k) is equal to tt/4 w 0.79. All other efficiencies are zeros. 

This natural cosine control law may however not match 
every desirable properties of the control law. For instance, if 
for some reason the control law has to be positive, one could 
chose 1 1 — y I ~\~ cos (|Aj — Afe|t) . The efficiency with respect 
to transition (J, k) is equal to 1/2. All other efficiencies are 
zeros. 

If there are no resonance, that is if Xj — Xk $ (A; — 
A m )Z for every {l,m} ^ {j,k}, one could consider the 
choice of a periodic Dirac pulse t i-> E„ e z <W|A.j-Afc|(*)> 
whose efficiency with respect to transition (j, k) is equal to 
one. If one desires to avoid the use of distributional control 
laws (Equation (|2]i has then to be understood in the measure 
sense), one may consider standard L 1 function close enough 



to Dirac pulses. Examples are presented in Section IV 



Finally, let us mention the algorithm described in [3] 
allows to build u* positive and piecewise constant. If Aj — 
Afe = 1 and the only resonances to be considered are such 



that |A; — A m | £ {(ii, a 2 , ■ ■ . , a p }, then the efficiency of u* 
with respect to transition (j, k) is 



n 



cos 



7T 

2a r . 



In the worst case, this algorithm guarantees E^ ,k \u) > 
Tln>2 cos(7r/2n) > 0.43, what may seem poor with respect 
to the cosine law. However, this algorithm is especially 
useful to handle the case of high order resonances. Indeed, 
if oi, CJ2, • ■ • , dp are all greater than N, then the efficiency 
with respect to transition (j, k) is greater than 



exp 



7T 

'In 



487V 3 



which tends to one as N tends to infinity. 

IV. Rotation of a planar molecule 

In this Section, we apply our results to the well studied 
example of the rotation of a planar molecule (see, for 
instance, [12], [4], [3]). 

A. Presentation of the model 

We consider a linear molecule with fixed length and center 
of mass. We assume that the molecule is constrained to stay 
in a fixed plane and that its only degree of freedom is the 
rotation, in the plane, around its center of mass. The state of 
the system at time t is described by a point 6 n- ip(t, 8) of 
L 2 (£l, C) where f2 = H/2nZ is the one dimensional torus. 
The Schrodinger equation writes 



i^(t, 0) = -A^(t, 0) + «(t) cos(0)V(i, 0), 



(6) 



where A is the Laplace-Beltrami operator on fi. The self- 
adjoint operator —A has purely discrete spectrum {k 2 , k £ 
N}. All its eigenvalues are double but zero which is simple. 
The eigenvalue zero is associated with the constant functions. 
The eigenvalue k 2 for k > is associated with the two 
eigenfunctions 9 H> cos(k8) and 8 H> ^=sin(A;0). The 
Hilbert space H — L 2 (fl, C) splits in two subspaces H e and 
H Q , the spaces of even and odd functions of H respectively. 
The spaces H e and H Q are stable under the dynamics of (I3|, 
hence no global controllability is to be expected in H. 

B. Non-resonant case 

We first focus on the space H a . The restriction A of iA 
to H Q is skew adjoint, with simple eigenvalues (— \k 2 )k^ 
associated to the eigenvectors 



0h-> 



1 



sin(fc0) 



fceN 



The restriction B of ip n- — icos(#)?/> to H a is skew-adjoint 
and bounded. The pair (A, B) satisfies Assumption [T] and is 
weakly-coupled (see [5, Section III.C]). 



The Galerkin approximations of A and B at order N are 
/ i ■•• \ 
4i '•. : 

: '■■ '■■ 

\ ••■ N 2 i J 

( 1/2 •• 
1/2 1/2 '■ 

o '■■ 



B (N) 



and 



\ 







V 



1/2 




1/2 






••• 1/2 

Our aim is to transfer the wave function from the first 
eigenspace to the second one. The numerical simulation 
will be done on some finite dimensional space C^. The 
controls we will use in the following have L 1 norm less 
than 13/3 and, from Proposition fo\ the \A\z norm of 
T"(0i) will remain less than exp(13/2) ps 665 for all 
time. From [5, Remark 4], the error made when replacing 
the original system by its Galerkin approximation of order 
^/l3e [ ^y3/^(F 2 288228 is smaller than e = 1(T 2 . This 
estimate is indeed very conservative and it can be improved 
using the regularity of the operator B. 

From [5, Section IV.C], for every integer I, for every 
t in [0, +oo), for every locally integrable control u (not 
necessarily periodic), 



< 



1 

fc! 



u{t)\At 

o 

13/3, then ||7r 22 i3(Id 



As a consequence, if \\uWl 1 < 
7i-22)TJ l (0i)|| < 5.10" 7 for every t in [0,+oo). Using this 
inequality, one gets that the error made when replacing the 
original system by its Galerkin approximation of order 22 is 
smaller than e = 3.10~ 6 when < 13/3. 

The transition between the levels 1 and 2 is resonant, 
indeed, 5 2 - 4 2 = 9 = 3(2 2 - l 2 ). Nevertheless, for 
every {hjh} ^ {1,2} such that — A/ 2 £ 3Z and 
(<f>l u B(f>i 2 ) 7^ 0, one has li > 2 and l 2 > 2. Hence, for 
every ^-periodic function u, the limit of the propagator 
X?£?(t, 0) leaves invariant the subspace generated by </>i 
and q>2 and the result of Theorem [5] applies (without having 
to check that all efficiencies of u for the transition (Zi,^) 
with h — l 2 £ 3(Z \ {1}) are zero). 

We illustrate the notion of efficiency on some examples 
of control, namely u* : t M> cos z (3t) for I £ {1, 2, 3, 4, 5}. 

The efficiency is zero when / is even. In numerical 
simulations, the quantity \{(j>2,X^ 22 ^{t,Q)(j)i)\ is less than 
2.1(T 5 for every t < 500 (see Figure 1 for I = 2). 

When I is odd, the efficiency is not zero. To estimate 
numerically the efficiency, one considers, for n £ {1, 10, 30}, 
the first maximum p 1 of 1 1— > \(<p2, -^nv) 0)^i)|> reached 
at time v, and computes 

(1 — p^nir 

2|(0 1 ,S0 2 )|/ o tt | U *(r)|dr' 



TABLE I 

Numerical Efficiencies of some periodic shapes 

1.0 



Control u* 


n 


Time ft 


Precision 


Numerical 


D.9 - 


(Efficiency) 






1-pt 


Efficiency 


0.8 - 




n = 1 


6.8 


2.10" 2 


73% 


0.7 - 


t i-> cos(3t) 


n = 10 


63 


4.10" 4 


78% 


0.6 - 


tt/4 w 79% 


n = 30 


189 


3.10" 5 


78% 


0.5 - 




n = 1 


8.9 


2.10" 2 


83% 


0.4- 


t m- cos(3t) 3 


n = 10 


84 


2.10 -4 


88% 


9yt/32 ft! 88% 


n = 30 


252 


2.10" 5 


88% 


0.3 - 




n = 1 


10 


7.10" 3 


93% 


0.2 - 


1 1 — ^ cos(3t) 5 


n = 10 


101 


2.10" 4 


92% 


o.i - 


75tt/256 ft 92% 


n = 30 


302 


2.10" 5 


92% 


o.o 





The Scilab source codes used for the simulation are 
available on the web page [7]. We sum up the results in 
Table 1. 



Fig. 2. Evolution of the square of the modulus of the second coordinate 
when applying the control t (—> cos 3 (3£)/30 on the planar molecule (odd 
subspace) with initial condition <j)\. The simulation has been done on a 
Galerkin approximation of size TV = 22. 
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Fig. 1. Evolution of the square of the modulus of the second coordinate 
when applying the control : t i— > cos 2 (3t)/30 on the planar molecule (odd 
subspace) with initial condition <j>\. The simulation has been done on a 
Galerkin approximation of size N = 22. 



C. Resonant case 

We focus on the space H e . The restriction A of iA to 
H e is skew adjoint, with simple eigenvalues (— ifc 2 ) 



fceNu{o} 
with (f)k \ i y 



o : a n- —0=. The restriction 

" Z7T 



associated to the eigenvectors 
--4= cos(fcf?) for k in N and 4 

B of ip i-> — icos(0)-0 to H e is skew-symmetric. The pair 
(A + i, B) satisfies Assumption [I] The translation from A to 
A + i induces just a phase shift and will be neglected in the 
following. 

The Galerkin approximation of A and B at order N are 



/ 
i 











and 



I 


l/s/2 







\ 




1/2 























■■• 1/2 





1/2 




V o 





1/2 





/ 



(N — 1) i / 



B (N) 



Our aim is to transfer the population from the first 
eigenspace, associated with eigenvalue 0, to the second one, 
associated with eigenvalue i. The transition (1, 2) is resonant 
(indeed 2 2 - l 2 = 3 = 3(1 2 - 2 )), and unlike what happens 
on the space of odd eigenfunctions, the limit matrix M' 
does not necessarily stabilize the space spanned by cf>i and 
<j>2 for every 27r-periodic function u*. Note however that B 
only connects level 2 to levels 1 and 3. In other words, it is 
enough to find a 27r-periodic function u* such that E^ 2 ' 3 \u*) 
is zero and E^^(u*) is not zero (and as large as possible) 
to induce the desired transfer. This is achieved, for instance, 
with the sequence of piecewise constant controls build in 
[3], for which the efficiency with respect to transition (1, 2) 
tends to cos(7r/6) and the efficiency with respect to transition 
(2, 3) is zero. Another example is presented on Figure 3. 

V. CONCLUSIONS AND FUTURE WORKS 

A. Conclusions 

The contribution of this paper is twofold. First, we have 
shown how simple regularity hypotheses can be used to 
approach with arbitrary precision an infinite dimensional 
system with its finite dimensional Galerkin approximations. 
Using this finite dimensional reduction, we then used clas- 
sical averaging techniques to obtain a proof of a well 
known experimental result about periodic control laws for the 
bilinear Schrodinger equation. As byproduct, we introduced 
the notion of efficiency, which characterizes the quality of 
the shape of a given control law. 




10 20 30 40 50 



Fig. 3. Evolution of the square of the modulus of the second coordinate 
when applying the control A cos(i) + jj on the planar molecule (even 
subspace) with initial condition <f>\. The simulation has been done on a 
Galerkin approximation of size N = 22. Precision 1 — p* is equal to 
2.10 -3 . Numerical efficiencies are 38% (theoretical: 3/8) for the transition 
(1,2) and less than 5.10 -4 for the transition (2,3) (theoretical: 0). 

B. Future Works 

Most of the points in this paper are merely a starting point 
to further investigations. Among other, we plan to study the 
generalization of the notion of weakly-coupled systems for 
systems with continuous or mixed spectrum. 
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